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Abstract 

The study of charmonium dissociation in heavy ion collisions is generally per- 
formed in the framework of effective Lagrangians with meson exchange. Some stud- 
ies are also developed with the intention of calculate form factors and coupling 
constants related with charmed and light mesons. These quantities are important 
in the evaluation of charmonium cross sections. In this paper we present a calcula- 
tion of the lvDD vertex that is a possible interaction vertex in some meson-exchange 
models spread in the literature. We used the standard method of QCD Sum Rules 
in order to obtain the vertex form factor as a function of the transferred momen- 
tum. Our results are compatible with the value of this vertex form factor (at zero 
momentum transfer) obtained in the vector-meson dominance model. 
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1 Introduction 



One of the most relevant topic in the physics of relativistic heavy ion collisions 
is the interaction of charmonium with nuclear matter, since the charmonium 
suppression is one of the most evident signals of formation of the quark-gluon 
plasma (QGP) [1]. The J/tp production and absorption in hadronic matter 
is still an open theoretical discussion [2]. It is important to understand the 
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two different ways of charmonium absorption: by the nucleons and by the 
co-mover light mesons (n, K, p,u, etc...). Calculations of form factors and 
coupling constants related with charmed and light mesons are of great impor- 
tance in evaluating charmonium cross sections [3] . Recent results of B AB AR, 
CLEO, BELLE and SELEX on D mesons spectroscopy stimulate calculations 
of physical quantities like vertices involving D mesons [4]. 

In the low energy regime (y/s < lOGeV), theoretical studies are performed 
mainly through the use of effective Lagrangians. In this context, one has a 
good control of the relevant symmetries that underly the dynamics of the 
process. On the other hand, it is necessary to know the values of form factors 
- associated with the vertices - with some precision. The choice of a lower 
or higher value of the form factor may change the final cross section in some 
orders of magnitude. 

In refs. [5-7] effective models for J/ip absorption in hadronic matter are pro- 
posed. The Lagrangians involve charmed and light pseudoscalar mesons (P) 
and also vector mesons (V). The interaction among these particles occurs 
through three-point (PPV and VVV) and four-point (PPVV and WW) 
vertices. One possible vertex in some of these models involves one u and two 
D mesons. On the other hand, in the framework of vector meson dominance 
(VMD) model, one can obtain the value of the coupling constant (at zero mo- 
mentum transfer) for such a vertex as g^DD ~ —2.84 (Appendix A of [7] and 
also [8]). 

The QCD sum rules (QCDSR) method [9] have been used in many works to 
calculate cross-sections, form factors and coupling constants (see for instance 
[10-14]). More recently, other works have used this method to obtain quantities 
related to charmonium suppression [15,16]. These works stimulate the use 
of QCDSR in these kind of problems. Some works have used form factors 
calculated with QCDSR method for charmed and light mesons [3] to study 
the J j ip dissociation. These form factors depend on the transferred momentum 
and this dependence is taken in to accout in the calculation of cross sections. 

Following these ideas, we present in this paper a study of ujDD vertex based 
on the QCDSR. This technique allows one to obtain hadronic quantities in 
terms of quark and gluon properties. Our results are compatible, within error 
estimates, with the value obtained through the VMD model. 

We follow the standard procedure of QCDSR. We calculate the Operator Prod- 
uct Expansion (OPE) and the phenomenological contributions for the corre- 
lation function of ooDD vertex and we equate both contributions following the 
principle of quark-hadron duality. In order to suppress higher order contribu- 
tions from the OPE side as well as higher resonances (and continuum) from 
the phenomenological side, we use the Borel transform in both sides of the 
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Fig. 1. The loDD vertex, with mesons momenta, in two situations: (a) D off-shell 
and (b) oj off-shell . 

equation, obtaining the sum rule. 

We perform the numerical integration of the sum rule to estimate the coupling 
constant. This coupling constant is a function not only of the transferred 
momentum Q 2 but also of the Borel masses. In general one considers the 
dependence of decay constants (fa and f u ) with Borel mass to improve the 
stability of the coupling constant with respect to the variation of the Borel 
masses. 



2 Sum Rules for the vertex 



Vertex with D off-shell. We start considering the ooDD vertex with one of 
the D mesons off-shell (fig. 1-a). The three-point correlator is given by 

r<?W,Q) = / d 4 xd 4 ye^e^yT^(x,y) , (1) 

where T^(x,y) = (0\T{j D (x)j' D (y)jj l (0)}\0); with the currents given by jo = 
iu^ 5 c for D and = u^^u for uj. 

In the QCD side we perform the operator product expansion (OPE) [17] of 
the correlator r( D ) (x, y) 



Tf\x,y) = A».l + B,.{qq) + ... (2) 

where 1 is the identity operator, A^(x,y) is the perturbative contribution, 
(qq) is the quark condensate and B^(x,y) is the respective coefficient. Next, 
we present the calculations of A^x, y) and B^x^y) with some detail. 



Apix, y) = Tr [Su^-y^S^y - x)~f 5 S u4 (x - 0)7 M ] , (3) 

where S Ut d is the light quark u (or d) propagator and S c is the quark c propaga- 
tor. After Fourier transforming the propagators, defining the variables: s = p 2 , 
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u = p' 2 , t = q 2 , using Cutkosky rules and some extensive calculation we get 
the dispersion relation 1 



Ay.{p,p',q) 



1 



CO CO 



ds I du- 



An 2 J J (s — p 2 )(u — p' 2 ) ' 
with the double discontinuity written as 
Vl D) =pp» + p'p',, 



(4) 



(5) 



and the spectral densities given by: 



P = 



3 (u — m 2 ) 6 



F, p' 



- —G , 



with 



^ 7r (s — t — u + 2m 2 ) 
F = -- — 

_ns(s-t-u + 2ml) 2 
~ 4 X 3 / 2 (s,t,u) 



1 - 



(s-t-u + 2m 2 c )(s -t + u) 
X(s, t, u) 



where A is the kinematical function: X(x, y, z) = (x — y — z) 2 — Ayz 



(6) 



(7) 



The integration limits of eq. (4) (s m j„ and u m i n ) are determined from the 
conditions: 



s > 0, u>t 
(s-t-u + 2m\) 2 < X(s, t, u) . 



(8) 



We go to the Euclidean space with the transformations: p 2 — > — P 2 , p' 2 — > — P' 2 
and q 2 = t — > — Q 2 . And after performing the double Borel transform in P 
and P', we obtain: 



CO CO 



B P 2 [B P ,2[A^\ = / rfs / 



e - S /M 2 e -«/M' 2 



(9) 



1 Notation: A^ = A^-l. 
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Now we consider the non-perturbative contributions. Here we shall concentrate 
only in the quark condensate contribution. The first contribution is 2 : 



B^{p,p',q) = J d 4 xd 4 ye^ x e-^B^(x,y) , 



where 



(10) 



B^(x,y) = Tr [- — (qq)) .J 5 S c (y - x) l5 S u , d (x)^ 



After some trivial manipulation we get 
p p z [q z — m z ) 



(12) 



It is not difficult to see that after going to the Euclidean space and taking the 
double Borel transform this contribution vanishes. 

The second contribution is 



Bj?\x,y) = Tr S Ujd (-y)^ 5 S c (y - x) ,^y- — {qq)jl„ 



that results in 



Bj?\p,p',q) = 



p 2 (p' 2 — m 2 ) 



{qq) ■ 



(13) 



(14) 



Going to the Euclidean space and taking the double Borel transform we get 



B 



p' 2 



Bp»[BV>\\ = -m c (qq)e 



-m 2 /M' 2 



(15) 



There is also a numerically negligible contribution from the charm condensate, 
which we will not take into account in this calculation. 

After the double Borel transformation the correlator (1), in the phenomeno- 
logical side, can be written as 



B 



B P ,.[Tf\p,p',Q)\ phe 



xe 



-ml/M 2 e -m 2 D /M' 2 ( _ 2 / + m D + m l + Q 2 

\ " " ">i 



* m 2 D + q 2 

pA+B P 2[B P n[}m\] ,(16) 



2 Notation: = B^ 
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where g^ D is the ooDD vertex, my^ and fy^ are, respectively, the mass and 
the decay constant of meson M (u or D). HR represents the higher masses 
(continuum) resonances, whose contribution will be discussed below. 

One important thing to be considered is the model to be adopted for the spec- 
tral density p^ R (s,t,u) function (HR term in the phenomenological side). We 
should remember that the interpolating fields couple not only with the funda- 
mental state but also with all particles with the same quantum numbers. We 
should also consider the contribution of those states in the phenomenological 
side. Admitting the quark-hadron duality, one assumes in general 

p^ R \s,t,u)^p^ PE \s,t,u)9(s-s )9(u-u ) , (17) 

where so and uo are the continuum thresholds. 

Proceeding this way, we have after the Borel transformation: 

oo oo 

B P , [Bpa [HR]] = - JL Jdsj dupf PE \s,t,u)e-^ M2 e-^ M ' 2 . (18) 

so uo 

In the case under consideration, it is not difficult to see that p ( -^ )PE \s,t,u) — > 

V^(s,t,u). Now we have all ingredients for computing of g^fji)(Q 2 ) 3 from 
the sum rule: (9) + (15) = (16). 

Vertex with cu off-shell. Now let us consider the case when the to meson is 
off-shell (fig. l-b). We follow a similar procedure of the previous section. The 
three-point correlator is written as 

T^(x,y) = (0\Tj D (x)jl(y)jU0)\0). (19) 

For the perturbative contribution we obtained the double discontinuity 
V^(s, u, t) = p(s,t,u).p^ + p'(s,t,u)p'^, with the the spectral densities: 



(s — u — t) 

p(s, t, u) = +3 ^ 3/2(s - j .[m 2 c (m 2 c - s - u + t) + su] (20) 

f/{s,t,u) = -3 ^ 3 ~^^ .[m 2 c (m 2 c -s-u + t) + su] . (21) 

The light quark condensate contributions (_B M ) vanishes after double Borel 
transformation. 



The superscript (D) indicates that D is off-shell. 
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In the phenomenological side, the correlator (after double Borel transforma- 
tion) is given by 



B 



p2 



9 { :Uq 2 ) 



Xe -ml/M* e -ml/M* + p J + Bp2 [£ P , 2 [HR]] . 



(22) 



Now it is straightforward to construct the sum rules and to calculate the value 
of vertex form factor g^DD- 



Decay constants and thresholds. In the numerical evaluation of sum rules 
one can take into account the dependence of the decay constants with the 
Borel mass. In general this procedure improves the stability of the vertex sum 
rules with respect to the increasing of Borel mass. From the sum rules for the 
correlators of the mesons we obtained the decay constants as functions of the 
Borel mass: 



f D (M 2 ) 



~M 2 
4^ 



;i_ e -VM 2 ) + 2 



m q (qq) 
M 2 



mz 



so 



m 



D 



8tt 2 



ds 



> ™ 2 cY p -s/M* _ inJ(l(j)( ,„ 



e m 2 D /M' 



(23) 
(24) 



As a common practice one usually assumes for the continuum threshold: s = 
(m r +A r ) 2 , where m r is the mass of the ground state resonance and A r ~ 0.5 — 
0.6GeV. For consistency we verified that the value of A r above is suitable, at 
least for u> and D mesons. Such verification can be done once f r (M 2 ) (actually 
f r (M 2 ,so)) stabilizes fast enough with growth of M 2 . Then we assume that 
f r {M 2 — > oo, so) = fr(so) tends to the phenomenological (experimental) value 
of f r . 

For instance, in the case of the uj meson we can observe from (23): 
,2/^2 \ s (m w + A w ) 2 

UM - oo, so) - j^y 2 = (27r)2 , (25) 



assuming ftf 3 *) « 0.20GeV and m w « 0.78GeV, we obtain A w « 0.48GeV. 

Repeating a similar calculation for eq. (24) with: fjj 0.17GeV, 
1.9GeV, m c w 1.2GeV and (qq) w (-0.23GeV) 3 , we estimate A D m 0.53GeV. 
We have seen that A r ps 0.6GeV led to better stability of sum rules for larger 
values of Borel mass. 
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p' structure 
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Fig. 2. Behavior of (— g^ D ) with M 2 for some fixed values of Q 2 . 



3 Numerical results and discussions 



As one can observe, in the sum rules above we have two structures: one related 
to Pn and other related to p' . For D off-shell the last structure is more stable 
(with respect to Borel mass) than the first one 4 . Therefore we concentrated 
our efforts in the analysis of p'^ structure. 

The first step of the numerical analysis is to consider the behavior of the form 
factor vertex with respect to the variations of Borel masses. One important 
point that appears in this step is how to handle with the Borel masses M 
and M'. It is more or less intuitive that the Borel mass must be of the same 
order of magnitude of the mass of the corresponding particle. Then, for D 
off-shell, we take the ratio M/M' ~ m^/mo, and for uj off-shell we have 
M/M' « m D /m D . 

In fig. 2 we show the curves corresponding to g^D as function of Borel mass for 
some fixed values of Q 2 . We can note that g^^jj stabilizes for M 2 > 15GeV 2 . A 
similar analysis is performed for D off-shell, and we noted that stabilizes 
for M 2 > lOGeV 2 . 

The second step is to fix the Borel mass in values for which the vertex form 
factor stabilizes, and study its dependence with Q 2 . In figure 3 we plot our 
results for both situations: D off-shell with M 2 = lOGeV 2 (□) and uj off-shell 
M 2 = 15GeV 2 (•). 

4 On the other hand, for u off-shell one can see that both structures are equal. 
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* o) off-shell -lvf= 15.0 Gel/ 2 

1 D off-shell -fif =10.0 GeV 2 
- Fit: -g'Z(tf) «31/(<f+11) 



^ 



D 



12 3 4 

d (GeV 1 ) 

Fig. 3. Behavior of {—g^DD{Q 2 )) for oj and D off-shell (p' structure). 
We fitted some functions to the curves of figure 3. For u off-shell we found: 
-^(Q 2 )-2.9e-« 2 / 3 - 2 , (26) 



and for D off-shell 



-9 { JS D (Q 2 ) 



31 



Q 2 + ll 



(27) 



From these fits we can estimate g^DD ~ —2.9 at zero momentum transfer, 
which perfectly agrees with the result of VMD model. 

In conclusion, we should say that the use of QCD Sum Rules approach to 
evaluate the g^DD vertex is completely acceptable and demostrates that it is 
still an open way to investigate other important hadronic quantities related 
with the formation of the quark-gluon plasma. 
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